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Qh The QQ semirelativistic interaction in QCD can be simply expressed in terms 



of the Wilson loop and its functional derivatives. In this approach we present 



the QQ potential up to order 1/m 2 using the expressions for the Wilson loop 



given by the Wilson Minimal Area Law (MAL), the Stochastic Vacuum Model 



(SVM) and Dual QCD (DQCD). We confirm the original results given in the 
different frameworks and obtain new contributions. In particular we calculate 
up to order 1/m 2 the complete velocity dependent potential in the SVM. This 
allows us to show that the MAL model is entirely contained in the SVM. 
We compare and discuss also the SVM and the DQCD potentials. It turns 
out that in these two very different models the spin-orbit potentials show up 
the same leading non-perturbative contributions and 1/r corrections in the 
long-range limit. 
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I. INTRODUCTION 



Since the pioneering paper of Wilson a real breakthrough opened in the treatment 
of quark states and in this framework a lot of work was devoted to the study of the heavy 
QQ. The challenge was understanding low energy QCD dynamics and hence confinement. 
The main characteristics of the heavy meson and baryon spectrum are simple and cleanly 
connected to expectation value of the QQ and 3Q potentials. The size of the b and c systems 
extends over distances where confinement already plays a relevant role (only toponium can be 
described purely in terms of one gluon exchange plus higher order perturbative corrections 
but, as well-known, we cannot access its spectrum); moreover, due to the mean value 
of the quark velocities, the leading relativistic corrections can be appreciated and usefully 
tested on the data. Furthermore, a good understanding of the heavy quark semirelativistic 
interaction is the first step towards relativistic generalization. 

At the static level, the linear confining QQ interaction, corresponding to a constant 
energy density (the string tension a) localized in a flux tube between the quarks, emerges in 
lattice formulation of QCD and is contained in all the existing confining models, e.g., Wilson 
area law, flux tube model and all kind of dielectric and dual models. This corresponds also to 
the static limit of the Buchmuller's picture of a rotating quark-antiquark state connected 
by a purely chromoelectric tube with a pure transverse velocity and with chromomagnetic 
field vanishing in the comoving system of the tube. In this picture it follows simply that the 
non-perturbative spin-interaction is given only by the Thomas precession term. 

The spin- dependent relativistic corrections were calculated first by Eichten, Feinberg 
and Gromes as a correction to the static limit (Wilson-Brown-Weisberger area law 
result). The potential is expressed in terms of average of electric and magnetic fields that 
can also be calculated on the lattice. The Eichten-Feinberg-Gromes results, at least in the 
long range behaviour, have been reproduced on the lattice H[7J (for a detailed discussion see 
Sec. 6). Recently the spin-dependent potential was also studied in the context of the Heavy 
Quark Effective Theory M. 
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In the literature relativistic generalizations of these results were attempted in a Bethe- 
Salpeter context by constructing a Bethe-Salpeter kernel which give back static and spin- 
dependent potentials. Using a simple convolution kernel (i.e. depending only on the mo- 
mentum transfer Q), this amounts to considering a Lorentz scalar proportional to 1/Q 4 . 
The velocity dependent relativistic corrections were also obtained but they are strongly de- 
pendent on the type of "instantaneous" approximation chosen to define the potential and 
on the gauge. These non-perturbative velocity dependent corrections destroy the agreement 
with the data |)-3 1] and give origin to the puzzle of how reconciling the spin-structure (i.e. 
the Lorentz nature of the kernel) with the velocity corrections in one Bethe-Salpeter kernel. 
In this paper we will not deal with this problem starting directly from the 1/m 2 expansion 
of the potential without any relativistic assumption. However a first step in its resolution 
seems to be the correct inclusion of the low energy dynamics also in the spin-independent 
1 / m 2 corrections. Moreover from the knowledge of these and the spin dependent corrections 
we will obtain some important insights on the nature of the kernel. 

Recently a method to obtain the complete 1/m 2 quark-antiquark (and 3 quarks) po- 
tential, based on the path integral representation of the Pauli-type quark propagator, was 
given in |12| (see also Hl3| , [T^| and [fL5Q. This formulation is gauge invariant. The potential 
is obtained as a function of a generalized Wilson loop (i.e. any kind of trajectory for the 
quark and the antiquark can appear) and its functional derivatives. These are all measur- 
able on the lattice. In short it was obtained a constituent quark semirelativistic interaction 
with coefficients determined by the non linear gluodynamics. This is the ideal framework 
in which to formulate hypothesis on the Wilson loop behaviour (and so on the confinement 
mechanism) to be checked on the lattice and on the experimental data. 

First, to evaluate the non-perturbative behaviour of the Wilson loop, a modified minimal 
area law (MAL) was used (see Sec. 3). This reproduces the Eichten-Feinberg-Gromes 
results jm and gives a velocity dependent potential proportional to the flux tube angular 
momentum squared, so that, by including velocity dependent corrections, a "string model" 



3 



emerges (see |TT| , [T6|| ). Also the velocity dependent potentials seem to agree with recent 
available lattice data [I7|. 

However, the MAL represents an extreme approximation that gives the correct result for 
very large interquark distances and does not give insight into open problems such as the 
relation between the non-perturbative structure of spin and velocity corrections. For these 
reasons we have taken into account two models of confinement, the stochastic vacuum model 
(SVM) and Dual QCD (DQCD) which both give an expression for the whole behaviour of the 
Wilson loop and contain the area law in the long distances limit. It is interesting to realize 
that both models reproduce essentially the perturbative plus MAL results respectively in the 
limit of short and long distances but produce also subleading corrections. These allows us 
to understand better the physical picture. For example in the case of the non-perturbative 
spin-orbit interaction it turns out that the magnetic term cancels in the area law limit (zero 
magnetic field in the comoving framework) but presents 1/r suppressed corrections in the 
other two models. 

A careful comparison between the SVM and DQCD corrections and an investigation of 
the approximations in which they coincide seem to be of great importance to the aim of 
understanding the low energy gluodynamics contained in the Wilson loop. 

The plan of the paper is the following one. In Sec. 2 we briefly sum up the definition of 
the semirelativistic potential and the notations. In Sec. 3 we collect the results obtained in 
the MAL model. In Sec. 4 we briefly present the SVM and use it to evaluate the potential 
in the context of Sec. 2. In particular we obtain also the SVM velocity dependent potential 
which is new. We show that it satisfies important identities and we give the short and long- 
range limits. In Sec. 5 we introduce the DQCD potential and discuss the long-range limit. 
In Sec. 6 we discuss our results in connection with the up to now available lattice data and 
draw some conclusions. 
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II. THE QUARK- ANTIQUARK POTENTIAL 



In |L2| a Foldy-Wouthuysen transformation on the quark-antiquark Green's function was 
done and the result was written as a Feynman path integral over particle and anti-particle 
coordinates and momenta of a Lagrangian depending only upon the spin, coordinates, and 
momenta of the quark and antiquark. Separating off the kinetic terms from this Lagrangian 
it was possible to identify the heavy quark potential Vqq (closed loops of light quark pairs 
and annihilation contributions were not included): 

£ dtV QQ = i \og(W(T)) - £ JL dx» [s\ ({F ¥ (x))) - ^S l / kr p k 3 ({F, r (x))) 

2 A„2 



(D»F ufl (x)))) - \ £ -^-T s f dx» [ dx'°S\S k r 



x ( ((F lfl (x)F ka (x'))) - ($„(*)» ((F ka (x>)))) , (2.i; 



where 



F^ = d,A u - d u A, + ig [A„, A v \ , = \e^F pry , (2.2) 

D v F ViX = d»F u , + ig[A", F vy \ , (2.3) 



W(T) = P exp 



ig j> dx^A^x) 



(2.4) 



and 



1 fVAe iSYM< ~ A) f(A) 

((f(m - JVAe^ A ^rP{f(A)exp[tgf r dx^A,(x)]} 

/ DA e iS ™W& P exp ^ dxM^i)] ' { } 

The closed loop Y is defined by the quark (anti-quark) trajectories zi(t) (z 2 (t)) running from 
yi to x x (x 2 to y 2 ) as t varies from the initial time £j to the final time tf. The quark (anti- 
quark) trajectories Zi(£) (z 2 (t)) define the world lines r x (r 2 ) running from £j to tf (tf to £j). 
The world lines Ti and T 2 , along with two straight-lines at fixed time connecting y x to y 2 
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and xi to x 2 , then make up the contour T (see Fig. 1). [] As usual A^x) = A a ^(x)\ a /2, Tr 
means the trace over color indices, P prescribes the ordering of the color matrices according 
to the direction fixed on the loop and Sym(A) is the Yang-Mills action including a gauge 
fixing term. 

As the 1/m 2 terms in Vqq are of two types, velocity dependent Vvd and spin dependent 
Vsd, we can identify in the full potential three type of contributions: 

Vqq = Vb + Vvd + V SD , (2.7) 

with Vo the static potential. 

The spin independent part of the potential, Vo + Vvd, is obtained in ( |2.1|) from the zero 
order and the quadratic terms in the expansion of log(PV(r)) for small velocities zq(t) = 



Pi/mi and z 2 (t) = p 2 /m 2 . In the notation of [|T|,|l9|] the terms arising from this expansion 
can be rearranged as: 



i\og(W{T)) = f f dt V (r(t)) + V VD (r(t)) , (2.8) 

V VD (r(t)) = — ( Pl • P2 V b (r) + fi Pl • p 2 - Pl r 2 P2 r ) V c (r)| 

■mim 2 I \3 r z J J Weyi 

+ E —2 \p 2 M r ) + (k 2 - Pj ' r 2 Vl ' r ) w} , (2.9) 

~, mi I J V3 J r 2 J J wevi 



j=l 3 



where r(i) = zi(t) — z 2 (t) and the symbol { }we y i stands for the Weyl ordering prescription 
among momentum and position variables [T3 . 



The spin dependent potential Vsd contains for each quark terms analogous to those ob- 
tained by making a Foldy-Wouthuysen transformation on the Dirac equation in an external 
field (where ((F^)) plays the role of the external field), along with an additional term Vss 
having the structure of a spin-spin interaction. We can then write 



1 As a consequence J r dx^f^x) = (— dt(fo(zj) — zj ■ f(zj)), where Zj = (t,Zj(t)). The 
factor ( — 1) J+ accounts for the fact that world line T 2 runs from tf to i;. We also use the notation 
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VsD — V LS + ^Thomas + ^Darwin + V<3S , (2-10) 

using a notation which indicates the physical significance of the individual terms (MAG 
denotes Magnetic). The correspondence between ([2.10 ) and ( [2.1| ) is given by 



/ f dtVg™ = dx*S$ ((F^x))) , (2.11) 

^Thomas = E A / dx» P ) ((F^x))) , (2.12) 
^Darwin = E J~2 l dx» {(D" F v ^{x))) , (2.13) 

Jti - =1 orrij JFj 

f t! dtV ss = -±E ^^ T s / ^ / ^ ( ((Fi IJ/ (x)F ka (x'))) 



-((FUx))) ((F ka (x')))) ■ (2.14) 



In the well-known Eichten and Feinberg notation |4j] and taking also into account the Darwin 
potential and similar contributions arising from the spin-spin interaction p~^JT~9f| , the terms 
in Vsd can be rearranged as 

i, »=K4 + ^) AW(r)+i, * (r)l 



H — (Li ■ S 2 — L 2 ■ Si) ~V 2 (r) + — ( S ' rS 2- - is, ■ S 2 ) V 3 (r) 

1 



-Si-S 2 y 4 (r), (2.15) 

with Lj = r x pj. It is not possible to identify directly each Eichten and Feinberg potential 
with the terms contained in eq. fl2.1| ) without making some assumptions on the Wilson loop. 
This will be the aim of the next sections. But some observations are just now possible. The 
contributions to A(V + V a ) come from Voarwin and from Vss with j = j'. In the case j ^ j', 
Vss contributes to the tensor term V3 and to the spin-spin term V4. Finally, V\ receives 
contributions from both the magnetic (V^ AG ) and the Thomas precession term (Vxhomas) 
while the contributions to V2 come only from the magnetic term. 
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Due to the Lorentz invariance properties of the Wilson loop some exact relations for the 
potentials Vi and V a , V e can be obtained. The first was given by Gromes || for the 
spin-related potentials 

^[V (r) + V 1 (r)-V 2 (r)] = 0, (2.16) 



and the other one by Barchielli, Brambilla and Prosperi [14] for the velocity-related potentials 



Vd(r) + ivi(r) + iv„(r) - = , (2.17) 

H(r) + iv«(r) + ^ = 0. (2.18) 

Since these relations are due to the Lorentz invariance they must be satisfied by any good 
choice of the Wilson loop approximated behaviour. 

Summarizing, the static and velocity dependent part of the potential are given in terms 
of the expansion of the Wilson loop average (W(T)), while the spin dependent potentials are 
given as a sum of terms depending upon the quark and antiquark spins, masses and momenta 
with coefficients which are expectation values of operators computed in presence of a moving 
quark- antiquark pair. These expectation values can be obtained as functional derivatives of 
\og(W(T)) with respect to the path, i.e. with respect to the quark trajectories Zi(t) or z 2 (t). 
In fact let us consider the change in (W(T)) induced by letting Zj(t) — > Zj(t) + 5zj(t) where 
6z?(t i )=6zf{t t ) = 0: 



Varying again the path 

g 2 (((F^F^z,)}} - ((F^UFM))) = -<^_* ((F^^))). (2.20) 

All contributions to the spin dependent part of the potential can be expressed as first and 
second variational derivatives of \og(W(T)). Therefore the whole quark-antiquark potential 



depends only on the assumed behaviour of (W(T)). In the next sections we will discuss 
some of these assumptions and give for each of them the explicit analytical expression of the 
potential. 



III. MINIMAL AREA LAW MODEL (MAL) 



In Ref. (T^JT^I CW(r)) was approximated by the sum of a perturbative part given at 
the leading order by the gluon propagator D^ v and a non-perturbative part given by the 
value of the minimal area of the deformed Wilson loop of fixed contour T plus a perimeter 
contribution V: 



i log(W(r)) = i log(iy(r)) bK + i \og(W(T)) 
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-~g 2 £ dx^j dx\ iD^Xi - x 2 ) + crS min + . 



(3-i; 



Denoting by u M = u M (s, t) the equation of any surface with contour V (s G [0,1], t 6 
[£i,£f], u°(s, t) = t, u(l,t) = zi(i), u(0,t) = z 2 (t) ) and defining ut = u - (u • n) n with 
n = (du/ds)\du/ds\~ 1 , we can write: 



Sri 



min \ dt [ ds 

Jti Jo 



min / dt [ ds 

Jti Jo 



du 



ds 



' du 11 du^ \ ( du^ du^ " 
dt dt J \ ds ds / 
1 

2^ 2 



+ 



' du^ du^ 
dt ds 



1 



du 

~dt 



T. 



(3.2) 



which coincides with the Nambu-Goto action. Up to the order 1/m 2 the minimal surface 
can be identified exactly (see App. B ref. 0]) with the surface spanned by the straight-line 
joining (t,zi(t)) to (t,Z2(t)) with ti < t < U. The generic point of this surface is 



""min = t U min = S Z X (t) + (1 - s) Z 2 (t) 



(3.3) 



with < s < 1 and zi(t) and z 2 (t) being the positions of the quark and the antiquark at 
the time t. Then, the exact expression for the minimal area at the order 1/m 2 in the MAL 
turns out to be 
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Smin = / dtr f ds [1 - (s Zit + (1 - s) Z 2T ) 2 ]2 

Jti JO 



(it r 



1 - - ( Z^ T + Z 2T + Zit ■ Z 2T ) 4 



(3.4) 



The perimeter term is given simply by 

V = |xi - x 2 | + |yi - y 2 | + YjJ dt^Jz^z 



n 1 1 



(3.5) 



and it is clear that we can neglect the time-independent perimeter contribution to the po- 
tential in the limit of big time interval tf —t- x . By expanding also eq. (|3.5| ) at the 1/m 2 order 
we have 



i\og{W{Y)) 



LR - 1 dt ar 



l- \ (z 2 T + z 2T + zit • Z2T 

D v 



c 



k 



l 



(3.6) 



For what concerns the perturbative part in the limit for large tf—t\ the only non- vanishing 
contribution to the Wilson loop is given by 

4 „ i*t rk 



«log(^(r)) SR = -~g 2 jT f dt x dt 2 z v 2 {t 2 ) iD^( Zl - z 2 ) . 



(3.7) 



In the infinite time limit this expression is still gauge invariant. Expanding z 2 (t 2 ) around 
ti it is possible to evaluate explicitly from eq. ( p.7|) the short-range potential up to a given 
order in the inverse of the mass. Self-energy terms are neglected. 

So, in this framework the following (MAL) static and velocity dependent potential were 
obtained: 



a. 



V = --- + ar + C 

3 r 



(3.8) 



and the explicit expressions for the potentials are: 



14 r = ar , 

y 1 9 r 9 

VJr) = --ar- -C , VJr 



Vc{r) = -- — 
3 r 

1 

— ar . 
6 



ar 



9 4 

These potentials fulfil the exact relations (|2.17 ) and ( |2.18|) 



(3.9) 
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Moreover by evaluating the functional derivatives for the Wilson loop, as given by eqs. 
( p. 19| )- ( pT20| ) , we obtain also the spin-dependent potentials 

AK(r) = 0, *. Vl {r) = -a, ±V 2 (r) = ^ , V 3 (r) = 4^ , V 4 (r) = fvro^r) . 

(3.10) 

These potentials reproduce the Eichten-Feinberg-Gromes results 0] and fulfil the Gromes 
relation ( |2.16| ). Notice that, as a consequence of the vanishing in this model of the long- 
range behaviour of the spin-spin potential Vss and the spin-orbit magnetic potential V^ AG , 
there is no long-range contribution to V2, V3 and V4. Instead V\ has only a non-perturbative 
long-range contribution, which comes from the Thomas precession potential (gig) . 

The MAL model strictly corresponds to the Buchmuller picture || where the magnetic 
field in the comoving system is taken to be equal to zero. Let us first notice that the perimeter 
contributions at the 1 / m 2 order can be simply absorbed in a redefinition of the quark masses 
rrij — > rrij + C/2 (for details see |14|])- Then let us consider the moving quark and antiquark 
connected by a chromoelectric flux tube and let us describe the flux tube as a string with 
pure transverse velocity v t . At the classical relativistic level the system is described by the 
flux tube Lagrangian PH-pTfl 

2 

C = -J2m J ^V]-a [ r 1 - vf , (3.11) 

j=l J o 

with v( = vit r'/r + V2t(l — r'/r). The semirelativistic limit of this Lagrangian gives back 
the non-perturbative part of the Vq and Vvd potential in the MAL model (notice that the 
minimal area law in the straight-line approximation is the configuration given by a straight 
flux tube) 0. The remarkable characteristics of the obtained VVd potential is the fact that it 
is proportional to the square of the angular momentum and so takes into account the energy 
and angular momentum of the string: 



"For a discussion of the relation between the two models in the path integral formulation see [12| 
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VV'd = ~TK~ • L 2 + L 2 • LO - £ ~tij . (3.12) 

Ylm\m2 r [ brrij r J 

Finally, the non-perturbative spin-dependent part of the potential in this intuitive flux tube 
picture simply comes from the Buchmiiller ansatz that the chromomagnetic field is zero in 
the comoving framework of the flux tube. 

We notice that even if V\ seems to arise from an effective Bethe-Salpeter kernel which is a 
scalar and depends only on the momentum transfer, a simple convolution kernel cannot repro- 
duce the correct velocity dependent potential Q3.12D or equivalently ( |3.9| ) p2| . Nevertheless 
the behaviour (|3.12|) seems to be important to reproduce the spectrum H-irT],|TB|,pB|,E3 . 



IV. STOCHASTIC VACUUM MODEL (SVM) 



The SVM (see |15|j25[| and for a review |26|]) in the context of heavy quark bound state 
gives a justification of the MAL model avoiding the artificial splitting of the Wilson loop in a 
perturbative and a non-perturbative part. It reproduces the flux tube distribution measured 
on the lattice |[27|. Moreover it allows to go beyond the MAL model in a systematic way 



(e.g. with the so-called perturbation theory in non-perturbative background [55]). The whole 
non-perturbative physics is factorized in some correlation function which can be calculated 
on the lattice. 

The starting point is to express the Wilson loop average (^(r)) via the non-Abelian 
Stokes theorem |29|,|3(| in terms of an integral over a surface S enclosed by the contour T, 



and then to perform a cluster expansion In order to allow lattice calculations all these 
quantities are given in the Euclidean metric. Some care must be payed in converting it in 
the Minkowskian metric before putting in eq. (|2.1|) . 



(W(T)) = (P exp [ig J^S^F^Xq] 
i=l J- J s 



dS HU:j {u j ){F l _ llUl {ui,XQ) . . . /•;,,,,,!>/;../•())), 



(4.1) 



(4.2) 
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The cumulants ( ) cum are defined in terms of average values over the gauge fields ( ): 

(F(l)) cum = (F(l)> , (F(l)F(2)) cum = (F(1)F(2)> - (F(1))<F(2)> , ... (4.3) 

and P F^ u (u, x ) = P exp ig f™ dx^A^ix) F pv {u) exp [ig f*° dx^A^x)} where x is an arbi- 
trary reference point on the surface S appearing in the non-Abelian Stokes theorem (|4.1|) . 
In general each cumulant depends on 5* and on xo, but, as the left-hand side of eq. (|4.1| ) 
does not, it is expected that in the full resummation of all the cumulants (right-hand side 
of eq. (f4.2|) ) this dependence will disappear [30|. To minimize the required cancellations S 



is chosen to be the minimal area surface. 

Equation ( |4.2| ) is exact. The first cumulant vanishes trivially. The second cumulant gives 
the first non-zero contribution to the cluster expansion ( |4.2| ). In the SVM one assumes that 
in the context of heavy quark bound states higher cumulants can be neglected and the second 
cumulant dominates the cluster expansion, or, in other words, that the vacuum fluctuations 
are of a Gaussian type: 

2 

log(W(r)} = - 9 - J dS^iv) J s dS Xp (v){F^{u,x )F Xp (v,x )) cum . (4.4) 

Neglecting the dependence on xq and on the arbitrary curves connecting xq with u and 
v which seems to be relegated to higher correlators, the Lorentz structure of the bilocal 
cumulant implies that it can be expressed as [|T^]: 

(F flu ( y U,X )F X p(v,Xo))cnm = (F^(u, x ) F Xp (v , x )) 

1 

+ 2 



d 

((u - v) x 5 up - (u - v) p 5 vX ) 



du p 



d 

+ dv~ " U ~ v ) pSflX ~( u ~ w )a<W 



D\[{u — v) ) (4.5) 



p_g 2 (TrF^(O)F^(O)) 



36 D(0) + D x (0) 

Eq. (f4.4|) and ( |4.5| ) define the SVM for heavy quarks. The correlator functions D and D\ are 
unknown. The perturbative part of Di, which is expected to be dominant in the short-range 
behaviour, can be obtained by means of the standard perturbation theory: 
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Df CTt (x 2 ) = ^ Qs — + higher orders. (4-6) 

Instead the only information which we know about the non-perturbative contributions to D 
and D\ come from lattice simulations. A good parametrization of the long-range behaviour 
of the bilocal correlators seems to be [[3^.[33|] : 



(3 D LR (x 2 ) = d e~ slxl , 6 = (1 ± 0.1) GeV , d = 0.073 GeV 4 , (4.7) 
P D\ K (x 2 ) = d 1 e~ 5l W , Si = (1 ± 0.1) GeV , d 1 = 0.0254 GeV 4 , (4.8) 

Up to order 1/m 2 the minimal area surface can be identified, as in the previous section, 
with the straight-line surface (|3.3|) . In particular, since dS^iu) = dt ds du^(t, s)/dt 
du v (t, s)/ds, we have: 

dSij(u) = dt ds Tj(t) ; 

dSij(u) = dt ds (s Zixit) + (1 — s) z i2 (t)) rj(t) . 

From flOD and (|4.5| ) and taking in account (|3.3| ) we have calculated explicitly \og{W(T)). 
Considering time interval much larger than the typical correlation length of D and D±, up 
to order 1/m 2 we have (for details see the appendix): 



V (r)=P 



J + °° dr jjT d\(r - A) D(t 2 + A 2 ) + jT dX^D^r 2 + A 2 )j , (4.9) 



r 



2 



2 D l {r 2 + r 2 ) \ , (4.11) 



P f +co , i r » / 2 3 A 1A2 4A3 3r2> 



/ dX (~X + D^r 2 + X 2 )) , (4.12) 



o 



4 4 r 4 r 
/ r A 2 2 A 
2 7-oo Jo 1 + V ~ 3 r 2 / 



K(r) = £ f °° rfr f rfA f- r - + - - |^ D(r 2 + A 2 ) . (4.13) 
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Result ( |4.9| ) was found in fI5 |, whereas fl4.10| )-( fLl3| ) are new. We note that these expressions 



for the potentials Vq and Vb, V e satisfy identically the Barchielli-Brambilla-Prosperi 
relations fl2.17| ) and fl2.18| ). Of particular interest seems to be the potential V e that has only 



non-perturbative contributions in the bilocal approximation. 

To evaluate the spin dependent part of the potential, the only terms which we need 
are those with one and two field strength insertions (taking in account that ((D u F up (x))) = 
d u ((F u ^x)))). By means of eq. ( gT| ) and (g2p and Q: 



/ + OO ( fT 1 I ~) 

^ drjjf d\-D{r 2 + \ 2 ) + -D l {r 2 + r 2 )^ 



g((Fu(zi))) 


= P 


{kin 


~ ZnTi 




+ P 




- z i2 ri 




= P 




- z i2 ri 




+ P 




- ZaTi 



+oo rr 1 / \ \ 

+ °° dr I f d\-,D(r* + A 2 ) + ilMr 2 + r 2 ) 
-oo [jo r z 2 

+oo rr 1 / \ \ 

dr / dX-(l--)D(T 2 + X 2 ) 

-oo jo r \ r / 



+oo 



o 



/(((F^(^)F Ap (z 2 ))) - ((F, u ( Zl )))((F Xp (z 2 

P(6^6 vp - S, p 5 uX ) (D(r 2 + r 2 ) + D^r 2 + 

d 

+P (r p r x 5vp - r^Tpbyx + r u r p 5 pX - r v r x 5 P p) -^D^t 2 + r 2 ) , 



r 4 = T = ti-t 2 



In this way we obtain the following expressions for the spin-dependent potentials in the SVM 
(confirming the results obtained in |l5j with a different derivation): 



AV a (r 



V 3 (r 
VAr 



self-energy terms , 



-P 



dr dX 1 - - D(t 2 + A 



P / + °° dr { f dX- D(r 2 + A 2 ) + Kd^t 2 + r 2 )) , 
J-oo [Jo r 2 J 



-P 



-oo Q 

*, dT r V 



L>i(r 2 + r 2 



y dr j 3D(r 2 + r 2 ) + SZ^r 2 + r 2 ) + 2r 2 — 73 



i(r 2 + r 2 ) 



(4.14) 
(4.15) 

(4.16) 
(4.17) 
(4.18) 
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Potentials ( [4.9| ), (|4.15[) and ( |4.16| ) satisfy identically the Gromes relation ( |2.16| ). An appli- 



cation of the spin potentials to the bb and cc spectrum, with a discussion on the different 
type of parametrization of the correlation functions, can be found in [53,^1]. 

In the short-range behaviour (r — > 0), assuming that all the relevant contributions come 
from the perturbative part of Di ( |4.6|) eqs. (|4.9| )- (|4.13| ) and (|4.15| )-( }4.14 ) exactly reproduce 



(after subtracting the self-energy contributions) the a s -depending part of eqs. (|3.8| ), (|3.9| ) 
and (|3.10| ) of the MAL model. We observe that no gauge choice is necessary in this approach, 
which is manifestly gauge- invariant. Moreover we note that the short-range behaviour of the 
Di correlator is not ad hoc but emerges straightforwardly from the comparison with the a s 
expansion of the Wilson loop. 

In the long-range behaviour (r — > oo): 

V {r) = a 2 r + ±cf - C 2 , ± Vl (r) = -a 2 + ^ , ±V 2 (r) = ^ , (4.19) 
2 dr r ar r 

V3 and V4 fall off exponentially and 
A\4(r) = self-energy terms , 



^2 



with 



/+00 /*oo 
dr / dX D{t 2 + A 2 ) , 
-00 Jo 

/~\~ OO p OO f' ~\~ 00 p 00 

dr d\\D(r 2 + \ 2 ), CP=P dr dX A DAt 2 + A 
-00 JO J-od JO 

/+00 /*oo 
dr / d\ A 2 D(t 2 + A 2 ) , 
-00 Jo 

/+00 roc 
dr / d\ A 3 D{t 2 + A 2 ) . 
-00 Jo 



2 \ , 



By means of parametrization ( |4.7| )- (|4~8| ) we have 

ivd Ad 
a 2 = — ~ 0.2 GeV 2 , C 2 = — ~ 0.3 GeV , 


^ = ^^0.1 GeV, D 2 = ^^0.7, E 2 = ~ 2.3 GeV" 1 , 
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and 



V 3 {r) = 1.2(6) dJ8 x ne~ Sl 



V 4 (r) = 1.2(6) -= r2e~ dr + 1.2(6) r5e" dir - 2.5(2) d x J 8 X r?e 

Identifying u% with a and C| /2 — Ci with C then at the leading order in r — > oo the 
spin-dependent and velocity-dependent SVM potentials reproduce the long-range behaviour 
of the potentials ( |3.8|) , ( [3.9|) and ( |3.10|) in the MAL model. Notice that the constant terms in 
the static and velocity dependent potentials turn out in the same combination as necessary 
to be reabsorbed in a redefinition of the quark masses. Some differences emerge at the 
next orders. In the SVM the magnetic contribution to the spin-orbit potential (which we 
called V^s AG in Sec. 2) is not exactly zero in the long-range behaviour but gives some 1/r 
corrections. For this reason the potential dV^jdr does not vanish and the potential dV\/dr 
presents a 1/r correction to the Thomas precession term. Notice, also, that in the SVM 
the tensor potential V3 and the spin-spin potential V4 are exponentially decreasing with the 
distance r but not identically zero as in the MAL model. In the next section we will see how 
the Dual QCD model is able to reproduce this behaviour. Finally a very rich structure of 
entirely non-perturbative 1/r and 1/r 2 corrections emerges in the velocity dependent part 
of the potential. A lattice study of this kind of contributions is in progress [lTj and in the 



light of eqs. (|4.20|) should give an interesting check on the validity of the stochastic vacuum 
approach in the velocity dependent sector of the potential and possibly some new indications 
on the behaviour of the correlator function D. A last comment on the fact that A14 is not 
r dependent. This is a direct consequence of the bilocal approximation which we have 
adopted. In principle nothing prevents us from the existence of r dependent contributions 
coming from higher order cumulants. We think it will be an important task to estimate such 
kind of contributions and compare it with lattice results (for a more detailed discussion see 
0). 
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V. DUAL QCD (DQCD) 



The duality assumption that the long distance physics of a Yang-Mills theory depending 
upon strong coupled gauge potentials is the same as the long distance physics of the 
dual theory describing the interactions of weakly coupled dual potentials = J2a=i C«A a /2 
and monopole fields <B< = ELi Bf\ a /2, forms the basis of DQCD M. The model is 



constructed as a concrete realization of the Mandelstam-t'Hooft |3B| dual superconductor 
mechanism of confinement. Indeed, the explicit form of the Lagrangian expressed in terms 
of the dual potentials is not known in a non-Abelian Yang-Mills theory. Since the main 
interest is solving such a theory in the long-distance regime, the Lagrangian £ e g- is explicitly 
constructed as the minimal dual gauge invariant extension of a quadratic Lagrangian with 
the further requisite to give a mass to the dual gluons (and to the monopole fields) via a 
spontaneous symmetry breaking of the dual gauge group. 

We denote by (W e f[{T)) the average over the fields of the Wilson loop of the dual theory 



JVC,VBVB 3 e*f d ^ GS ^ +c ^ 

where £gf is a gauge fixing term and the effective dual Lagrangian in presence of quarks is 
given by 

C cS (Gl) = 2 Tr {-\g^G„ + IW) 2 } ~ U(BA • (5-2) 
U{Bi) is the Higgs potential with a minimum at a non-zero value £> i — B X 7 , B 02 = —B X 5 



and £>03 — B \ 2 . It was also taken B\ = B 2 = B. In ( |5.1| ) we have taken the dual potential 



3 The name Dual QCD has historical reasons, but can give rise to some confusions. We emphasize 
that the duality assumption concern only the long distance physics of a strongly coupled Yang-Mills 
theory as the gluonic sector of QCD. 
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proportional to the hypercharge matrix C M = C^Y Moreover 

V^Bi = d^Bi + ie[C^ Bt], e = Y' (5 " 3) 
g v — ff) n — ft n _u \ v 



dpC v - d u C, + G%,J y , (5.4) 
G %( x ) = gz^uap jdsj dr——5(x - y(s, r)) , (5.5) 

and ?/(s, r) is a world sheet with boundary T swept out by the Dirac string. Notice that 
dual potentials couple to electric color charge like ordinary potentials couple to monopoles 

The functional integral (W e fr(r)) determines in DQCD the same physical quantity as 
(W(T)) in QCD. The coupling in C e s{G^ v ) of the dual potentials to the Dirac string plays 
the role in the expression Q of the Wilson loop W{Y) of QCD Q in (W(r)). The 
assumption that the dual theory describes the long distance QQ interaction in QCD then 
takes the form: 

(W(r)) = (W eff (r)) , for large loops T. (5.6) 

Large loop means that the size R of the loop is large compared to the inverse mass (M _1 ~ 
(600 MeV) -1 ) of the Higgs particle (monopole field). Furthermore, since the dual theory is 
weakly coupled at large distances, we can evaluate (W c s(T)) via a semiclassical expansion 
to which the classical configuration of dual potentials and monopoles gives the leading con- 
tribution. This then allows us to picture heavy quarks (or constituent quarks) as sources 
of a long distance classical field of dual gluons determining the heavy quark potential. We 



mention here that DQCD reproduces the lattice flux tube distribution [38]. 

Eq. ( |5.6| ) defines the DQCD model for heavy quark bound states. Replacing (W(T)) by 
(W / efr(r)) in eq. (2.1) we obtain expressions for Vq and Vvd and by considering the variation 



'Doing so £ fj without quark sources generates classical equations of motion with solutions dual 



to the Abrikosov-Nielsen-Olesen magnetic vortex solutions in a superconductor [18,19]. 
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in (Weff(r)) produced by the change G^(x) — > G s „ v (x) + 5G^ u (x) we obtain also the field 

averages in terms of dual quantities: 

nUF ( 7 )\\-( lV -+i ^ M^ eff (H) 4 . F SihgjWesiT)) 

g{{F^( Zj ))) - (-1) - -g((G^( Zl ))) cS - (-1) -e^ Xa . 

(5.7) 

This gives a correspondence between local quantities in the Yang-Mills theory and in the 
dual theory. A similar expression can be obtained for the double field strength insertion in 

The weak coupling of the dual theory permits the explicit evaluation of (W e s{T)) by 
means of the classical approximation. Hence we have 

% log<W eff (r)) = - j dx Crt{G%) , (5.8) 

with C e s(G^ v ) evaluated at the solution of the classical equations of motion: 

d a (d a Cp - d p C a ) = -d a G% + jM° N , (5.9) 
( d , + ieC»fB = -~^, (5.10) 

where j^ ION = —6 e 2 C IM B 2 is the monopole current. The Dirac string is chosen to be a 
straight line connecting Q and Q since this is the configuration having the minimum field 
energy. As a consequence of the classical approximation all quantities in brackets are replaced 
by their classical values ((G>j,(:r)))efr = G flv (x) which are obtained by solving numerically the 
non- linear equations (|5.9|) - (|5.11|) . An interpolation of the numerical results for the potentials 
can be found in [Tj| (in particular in the first of these references it is possible to find also an 
application of the DQCD potentials to the heavy quarkonia spectrum). In the following we 
will give and discuss only the large distances limit of these potentials. 

In the long-range behaviour (r — > oo) the interpolation of Ref. []18| gives 

V Q {r) = or - 0.646v^ , (5.12) 
20 



d Jr/ , 0.681 , . . 

V x {r) = -a+ yfra s , 5.13 

ar r 

d jr . . 0.681 / ft . 1A . 

-rV 2 {r) = y/Wa B , 5.14 

ar r 



and 



V h (r) = -0.097 ar - 0.226^/^ , (5.15) 

V c (r) = -0.U6 ar -0.516 ^a~ s , (5.16) 

V d (r) = -0.118 or + 0.275^^, (5.17) 

7 e (r) = -0.177 ar + 0.258v^a;. (5.18) 

For the spin-spin interaction and for large distances it is possible to give the exact analytical 
expression of the potentials: 

^ 3 (r) = -a s M 2 + -M+4 , (5.19) 

4 p-Mr 

V 4 (r) = -a B M 2 . (5.20) 

3 r 

While V a is, at the moment, lacking either in an analytical or a numerical evaluation, and is 
formally given by ||19|| : 

A 3 jl 

(5.21) 

X=x' = Zi 



AK(r) = -Ay-(r) - E X^77^ p (x, x') 



where the first term is the color electric contribution to V a (V r NP (r) is the non-perturbative 
part of the static potential, so that V & is determined by the non-perturbative gluodynamics) 
and the second is the color magnetic contribution. G NP satisfies the equation 

(-A + 6 e 2 B 2 )G NP = - 6e ; 2jg2(x) , . (5.22) 

47r|x — x'| 

The potentials depend on the two free parameters a s = n/e 2 and a. In [|T8| the values 

(x = 0.18GeV 2 , a s = 0.39, (5.23) 

were used. The dual gluon mass M is related to this two parameters and is approximately 
given by: 
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M 2 ~ j— ~ (600 MeV) 2 . (5.24) 



Finally we observe that all these potentials satisfy identically the Gromes relation Q2.16| ) 
and the equivalent relations for the velocity dependent potentials ( |2.17| ) and (|2.18|) . 

From the comparison of eqs. ( |5.12|) - ([5.13|) with ( |4.19| ) if follows immediately that in the 
long-range behaviour the static and the spin-orbit potentials coincide completely in DQCD 
and in the SVM. Very important seems to be the agreement, which we note here for the 
first time, between the 1/r corrections in the two models. These corrections come from 
the physics beyond the minimal area law assumption and in fact there are not present in 
the MAL model (see (|3.10|) ). The coefficient of the 1/r contribution in dV\/dr and c/V^/dr 
is the same in DQCD and SVM and in both cases compatible with the constant term in 
the static potential Vq. The little difference between the constant in dVi/dr, dV^/dr and 
Vo can be understood in the SVM language as due to the presence of the small positive 
constant C% /2. The spin-spin interaction falls off exponentially in both the models. In 
DQCD the behaviour is like a Yukawa interaction, while eqs. ( |4.17| )-( |4.18| ) seem not to 



reproduce this behaviour at least with parametrization ( |4.7| )-( |4lj| ). This is, at the moment, 
an important disagreement because one of the basic feature of DQCD is that the magnetic 
interaction (like in the spin-spin case) is carried by a massive particle. Differences arise for 
large distances also in the velocity dependent sector and with respect to the MAL model. 
The factors in front of the or leading contributions to \4, V e are slightly different from 



those of eqs. (|3.9| ). The potentials Vt>, V c and V e present some additional constant terms 
which do not arise from the area law. Finally there are not 1/r corrections as in the SVM. 
Some of these discrepancies can be interpreted as due to a finite thickness of the flux tube 



in DQCD opposite to the infinitely thin flux tube in the MAL model |L9|]. Therefore in the 
two models the flux tube will have a different moment of inertia and give slightly different 
contributions to the velocity dependent potential. It is possible that these discrepancies 
will disappear if including higher order cumulants contributions in the SVM predictions. 
Other differences between the predictions of the two methods could have origin from the 
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very delicate interpolating procedure of the numerical solutions of the DQCD non-linear 
equations. The very soon available lattice results on the velocity dependent potentials JH] 
will possibly clarify the situation. 

VI. DISCUSSION AND CONCLUSIONS 

Using the same gauge invariant and physically transparent approach to calculate the 
complete semirelativistic quark-antiquark interaction for three different models (MAL, SVM 
and DQCD) we have shown the following points. 

• We have obtained the velocity dependent corrections in the SVM model which are new 
and present an interesting non-perturbative structure. 

• We have demonstrated that the minimal area law model is exactly reproduced in both 
the spin dependent and the velocity dependent sector of the potential by the long-range 
behaviour of the stochastic vacuum model. From now we can consider the MAL model 
simply as the r — > oo limit of the SVM for heavy quarks. Moreover this limit realizes 
also the intuitive Buchmiiller's picture of zero magnetic field in the flux tube comoving 
system. 

• In the spin dependent sector of the potential, both the SVM and DQCD not only 
reproduce the long-range behaviour given by the area law, but also give 1/r corrections 
to dVi/dr and dV%jdr. These corrections are equal in both models and very near to 
the absolute value of the constant term in the static potential (the SVM also supplies 
for the explication of this fact). This perfect agreement is absolutely not trivial and 
seems to be very meaningful, since it arises from two very different models in a region 
of distances in which the physics cannot be described by the area law alone. This 
is also remarkable to understand the kind of effective kernel that would describe the 
non-perturbative bound-states of constituent quarks. For example, it seems now clear 
that the vanishing of the magnetic part, given by the field average of eq. (|2.11|) , in the 
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non-perturbative region takes place only at the leading level in the long-range limit. 
Therefore, working in a Bethe-Salpeter context, there is no need to assume an effective 
pure convolution kernel which is a Lorentz scalar (a recent proposed Bethe-Salpeter 
kernel can be found in 



Velocity dependent contributions to the quark- antiquark potential are important. In 
fact the string behaviour of the non-perturbative interaction shows up when we consider 
the velocity dependent part of the potential p3| , |i"U| and this is also what the data require 
[p3fl . The derivation of the velocity dependent part using equation ( [2~1~D and the S VM 
is completely gauge invariant and seems not to suffer from the problems connected 
with the strong reduction dependence of the potentials obtained from Bethe-Salpeter 
kernels. In this way we reproduce the area law results and give a lot of new 1/r 
and 1/r 2 corrections, suppressed in the long-range behaviour. The velocity dependent 
structure which arises from the DQCD model differs slightly in the coefficients with 
respect to the area law behaviour. The main reason seems to be that the flux tube in 
DQCD has a finite thickness. It is possible that higher order cumulants can reabsorb 
this difference. 

The spin dependent potentials have first been evaluated on the lattice. The data in 
M confirm the long-range behaviour given in ( |3.10| ) and contained also in ( |4. 19| ) and 



( |5.12| )-( |5.14| ). Recent data show up the same long-range behaviour and do not yet 



allow to distinguish between parametrizations which differ at the next-to-leading order 
in the distance r. However they contain more information about the short-range region 
of the interaction (typically below the correlation length of 0.2 fm). Generally the data 
reproduce the perturbative results (which at the first order in a s can be read from ( [3.10|) 



putting a equal 0). The only exception is given by the short distance behaviour of 
dVi/dr which seems to be negative and proportional to 1/ r 2 . This contradicts the order 
a 2 calculation of the QQ potential (which contains the first non vanishing perturbative 
contribution to dV\/dr) given for example by Pantaleone and Tye [^0|. The reason 
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of this discrepancy could be explained by higher order perturbative contributions or 
by some at the moment unknown short range non-perturbative contribution (in the 
language of the SVM this contribution could arise from the correlation function D; 
an investigation in this sense of the recent short-range data on D given in the last 
reference quoted in |32j is going on). The problem is still open. Only recently some 



data on the velocity dependent potentials appeared |L7y32[- Probably more accurate 
data will be available in the next months. These results seem to confirm the long-range 
behaviour contained in (|3.9| ) (a dependent terms). More interesting is the case of the 
potential AV a which appears to be different from zero for r — > oo and show up a 1/r 
short-range behaviour. This behaviour has been recently explained in terms of SVM 
and DQCD @. 

In conclusion SVM and DQCD reproduce the flux tube distribution measured on the 
lattice and the general features coming from the area law. Both give analytical expressions 
for the Wilson loop (eqs. Q4.4D and ( |5.1| )) which describe the evolving behaviour of (W(T)) 
from the short to the long distances (we note that this can be useful in many different 
applications, see e.g. [fy]]) and both give some predictions which go beyond the asymptotic 
behaviour. But not all predictions are equal in the two models in the intermediate distances 
region, in particular in the velocity dependent sector of the potential, but also in the spin- 
spin interaction. Therefore, new lattice data sensitive to such kind of corrections seem to be 
urgent. Finally, work is in progress in evaluating the correlation function D and D\ in the 
DQCD context and in producing an extensive phenomenological analysis of the contribution 
of the new obtained potentials to the heavy and heavy-light quark spectrum. 
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APPENDIX: 

In this appendix we derive the static potential in the SVM (eq. (|4.9| )). The same 
technique was used to obtain the other potentials. Since the velocity dependent potentials 
involve long and tedious calculations, a program of symbolic manipulations was used in that 
case 

From eqs. ( }4.4| ), ( |4.5| ) and in the straight-line parametrization of the surface, it follows 
that: 
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log(W(T)) = ~ [ dS 4l (u) [ dS Aj (v) 
2 Js Js 



5^ (D (t 2 + (u — v) 2 + D 1 r 2 + (u - v 



d 



+(r dij + {u- v)i{u - v) j )—D l [t z + (u- v 



dr 2 ' 



+ 0(z(,z 2 2 ), (Al) 



with dS^iiu) = dtidsiri(ti), dS^iv) = dt2ds2Tj{t2) and r = t\ —t%. Expanding the functions 
of t 2 around t\\ 



r i(*a) =r j {t l )-r j {t l )T+ ... , 
(u - v)j = z 2j (t 1 ) - z 2 j(t 2 ) + sirj(ti) - s 2 r j (t2) = (s x - s 2 )r j (t 1 ) + ... , 

and taking for simplicity rj(t\) = Tj and A = s 2 — Si, we obtain 

(3 ftf rt f rJ rJ 



log(W{T)) = / ' dh ! ' dt 2 f dai f ds 2 r irj S i5 (d(t 2 + A 2 r 2 ) + D 1 (r 2 + A 2 
2 Jti Jti Jo Jo K 

+ 0{zlzl). (A2) 



d 



+ (r 2 5 lJ + X\r 3 )—D l (r 2 + X 2 r 2 ) 



Since: 



dsi [ ds 2 f ((s 2 - 
o Jo v 



1 rl-si rl 

2\ 



ds x / dX f(X ) = 2 dX (1 - A) f(X 2 ) 

J-si Jo 



we can write 
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log(Mr)) = -p [ k dh t dt 2 f 1 dX (1 - A) r 2 
Jti Jti Jo 



d 

'dT 2 ' 



+ ( r 2 + A 2 r2) - jDi(r 2 + A 2 r 2 



D(r 2 + X 2 r 2 ) + D^r 2 + X 2 r 2 ) 

+ 0(zl,z 2 2 ). (A3) 



Replacing rX — > A and taking in account that the time variables in (|A3|) are in an Euclidian 
space while the equation for the potential Q2.1D is in Minkowski, the static potential is given 
by 



+00 rr 

V {r) = (3 I dr dX (r - X) D{r 2 + A 2 ) 
-00 Jo 



/+OO fT 
dr / dX{r- A) 
-00 Jo 



d 



D 1 (r 2 + A 2 ) + (r 2 + A 2 )^ 1 (r 2 + A 2 ) 



(A4) 



where, also, the large time limit was performed. Finally, the identities 



dr r'-^D^r 2 + X 2 ) = U dr r^D x [r 2 + A 2 ) = -- / dr D^r 2 + A 2 ^ 

-00 ST Z J— 00 CLT Z J— 00 



and 



f dX (r - A) X 2 ^-Mr 2 + A 2 ) = \ f dX (r - A) X^-D^ 
Jo dr z Z Jo dX 



A 2 ) 



_I f dA ( r _2A) -Di(r 2 + A 2 ) , 
Z Jo 



give back the static potential in the form of eq. (|4.9| ). Taking in account the 0(zf, z 2 ) contri- 
butions in ( |A1|) and in the following equations, we obtain the velocity dependent potentials 
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